Traits and Characteristics of Interacting Dirac fermions 
in Monolayer and Bilayer Graphene* 
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The relativistic-like behavior of electrons in graphene significantly influences the interaction prop- 
erties of these electrons in a quantizing magnetic fleld, resulting in more stable fractional quantum 
Hall effect states as compared to those in conventional (non-relativistic) semiconductor systems. In 
bilayer graphene the interaction strength can be controlled by a bias voltage and by the orientation 
of the magnetic field. The flnite bias voltage between the graphene monolayers can in fact, enhance 
the interaction strength in a given Landau level. As a function of the bias voltage, a graphene 
bilayer system shows transitions from a state with weak electron-electron interactions to a state 
with strong interactions. Interestingly, the in-plane component of a tilted magnetic fleld can also 
alter the interaction strength in bilayer graphene. We also discuss the nature of the Pfafflan state 
in bilayer graphene and demonstrate that the stability of this state can be greatly enhanced by 
applying an in-plane magnetic fleld. 
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I. INTRODUCTION 

Graphene is a monolayer of carbon atoms, which has a 
two-dimensional (2D) honeycomb crystal structure. The 
unique feature of graphene is that the single-electron low- 
energy dispersion has the relativistic massless form com- 
monly attributed to the Dirac fermions, and the corre- 
sponding electron wave functions have a chiral nature 
The electronic band structure, first derived by 
Wallace in 1947 has two valleys at two inequivalent 
corners, K = i2TT/a){\,^) and K' = (27r/a)(|, 0), of 
the Brillouin zone, where a = 0.246 nm is the lattice 
constant. The low-energy dispersion at the valleys ^ = 1 
(/C- valley) and ^ = — 1 (X'-valley) is determined by the 
following relativistic massless Hamiltonian 0-S S Q 



p_ 
P+ 



(1) 



where p_ — — iPy, Pj^ = Pa; + i-Py; ^''^d p is the two- 
dimensional momentum of an electron. Here v-p sa 10^ 
m/s is the Fermi velocity, which is related to the hopping 
integral between the nearest neighbor sites. The honey- 
comb lattice of graphene consists of two sublattices A and 
B and the two component wave functions corresponding 
to the Hamiltonian ([T]) can be expressed as (i/'a; V'b)"^ ^o^' 
valley K and ('i/'s, ^a)^ ^^r valley K' , where -0^ and ip^ 
are wave functions of sublattices A and B, respectively. 

The two components of the wave function correspond 
to the quantum mechanical amplitudes of finding the 
'Dirac fermion' on one of the two sublattices. This sub- 
lattice degree of freedom is often referred to as pseudospin 



of Dirac fermions in graphene. It is directed along the di- 
rection of motion of the Dirac fermion in the conduction 
band, and opposite to the motion in the valence band. 
In other words, particles in graphene have opposite chi- 
rality in the K and K' valleys. The electron and hole 
wave functions are eigenfunctions of the helicity (chiral- 
ity) operator Physically, a certain direction of the 
pseudospin in the graphene plane corresponds to a ro- 
tation of the relative phases of the two components of 
the spinor eigenstates along that direction of motion of 
the Dirac fermion. The sublattice pseudospin chirality of 
Dirac fermions does not allow perfect backscattering (be- 
tween states of opposite momentum and opposite pseu- 
dospin) that has important consequences on the physical 
characteristics of graphene [1, 0] . 

The single-electron states obtained from the Hamilto- 
nian ([T]) has a linear relativistic dispersion relation of the 
form 



e(p) = ±t)pp, 



(2) 
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where the signs and '-' correspond to the conduction 
and valence bands, respectively. Each energy level ([2]) 
is four-fold degenerate due to two-fold spin and two-fold 
valley degeneracies. 



II. DIRAC FERMIONS IN MAGNETIC FIELDS 

In a magnetic field applied perpendicular to the 
graphene plane, the relativistic energy dispersion rela- 
tion ([2]) brings in very specific form of Landau levels of 
electrons in graphene [10, 11], and as a consequence the 
2D system displays unconventional quantum Hall effects. 
The Landau levels of electrons in graphene can be found 
from the Hamiltonian ([T|) by replacing the electron mo- 
mentum p with the generalized momentum tt = p + eA/c. 
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Here A is the vector potential. Then the Hamiltonian 
of an electron in a magnetic field perpendicular to the 
graphene monolayer in valley ^ takes the form 



(3) 



The eigenfunctions of the Hamiltonian ([3]) can be ex- 
pressed in terms of the conventional Landau level wave 
functions, (/)„ ^, for a particle obeying the parabolic dis- 
persion relation with the Landau index n and intra- 
Landau index m, which depends on the choice of the 
gauge. For example, in the Landau gauge {A^ — and 
Ay = Bx) the index m is the y-component of the mo- 
mentum, while in the symmetric gauge {A — x f) the 
index m is the z-component of electron angular momen- 
tum. For these wave functions, ^, the operators tt_^_ 
and 7r_ are the raising and lowering operators, respec- 
tively. 

The Landau eigenfunctions of the Hamiltonian (|3]) are 
then written in the form 



K 



i'B 



sgn(n)il"l ^0, 



(4) 



for valley K {^^ = \) and 



n.m 



sgn(n)il"l 

il"l(/), , 

^ I n I , T?i 



for valley K' = -1). Here C„ 
C„ = 1/'n/2 for n 7^ and 




sgn(n) 



(n = 0) 

1 > 0) , 
-1 {n < 0). 



(6) 



where positive and negative values of n correspond to the 
conduction and valence bands, respectively. The corre- 
sponding Landau energy spectrum takes the form 



(7) 



where ujg = \/2wp/^Q and Iq = y^h/eB is the magnetic 
length. 

The specific feature of the Landau levels in graphene 
is their square-root dependence on both the magnetic 
field B and the Landau level index n. This behavior 
is different from that in conventional (non-relativistic) 
semiconductor 2D system with the parabolic dispersion 
relation, for which the energy spectrum has a linear de- 
pendence on both the magnetic field and the Landau level 
index, e„ — fujjg{n -|- 1/2). In Fig. [T]the Landau levels 
in graphene are shown as a function of the perpendicular 
magnetic field, where the positive and negative Landau 
level indices (n) correspond to the conduction and va- 
lence bands, respectively. The Dirac nature of the elec- 
tron dynamics in graphene and the unique behavior of 




10 20 30 

Magnetic field (Tesia) 



FIG. 1: The Landau levels as a function of the perpendicu- 
lar magnetic field for a graphene monolayer. Numbers next 
to the lines are the Landau level indices. The Landau levels 
with positive energies (positive index n) and negative ener- 
gies (negative index n) correspond to the conduction and va- 
lence bands, respectively. The fractional quantum Hall effect 
(FQHE) [19[ that is discussed below, can be observed only in 
the Landau levels shown by the red and blue lines (n = ±1 
and 0). The strongest electron-electron interactions and cor- 
respondingly the more stable FQHE states are realized in the 
Landau levels shown by red lines. 



Landau levels in a graphene monolayer, were experimen- 
tally confirmed by observation of quantum Hall plateaus 
at fining factors v = A{n+\) [H, 13. 

The interaction properties of electrons within a single 
Landau level, i.e., disregarding the mixture of Landau 
levels due to electron-electron interactions, are entirely 
determined by the pseudopotentials Vm proposed by 
Haldane [3, which are defined as the energy of two 
electrons with relative angular momentum m. They are 
determined by the structure of the wave functions of the 
corresponding Landau level and for the n-th Landau level 
can be evaluated from the following expression 



^qV{q) [F^{q)fL^{q^ 



(8) 



where L^{x) are the Laguerre polynomials, V{q) = 
2716^ / {nlgq) is the Coulomb interaction in the momen- 
tum space, K is the dielectric constant, and F^^{q) is the 
form factor of the n-th Landau level. In what follows, 
all pseudopotentials are given in units of the Coulomb 
energy, Sfj = b^/k^q. Equation ^ is valid for all types 
of electron systems (non-relativistic, monolayer and bi- 
layer graphene, etc.) with well defined 2D Landau levels. 
The difference between these systems is in the expression 
of the form factors, F^{q). In a non-relativistic system, 
for which the Landau level wave functions are the 
form factors are obtained from F„((?) = i„ (q^/2) . In the 
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case of graphene, the n-th Landau level wave functions 
are given by E qs. ( IH-dSl), which results in the following 
expressions [TtI. I18( for the corresponding form factors 



(9) 

(10) 
for 



TABLE I: Characteristics of the pseudopotentials 
graphene and for conventional electron systems 



for 



K-i ( y 



With these form factors the pseudopotentials 
graphene are then evaluated from Eq. 

One of the unique manifestations of electron-electron 
interactions within a single Landau level is the formation 
of incompressible fractional quantum Hall effect (FQHE) 
states [l9[, which are characterized by a finite excita- 
tion gap, determined by the electron-electron interac- 
tions. These states are realized at the fractional filling 
of a given Landau level, e.g., at filling factors v = 1/m, 
where m is an odd integer. The properties of the FQHE 
states are completely determined by the corresponding 

(n) 

pseudopotentials Vm ■ The stability of the incompress- 
ible FQHE state, i.e., the magnitude of the FQHE gap 
depends on how fast the pseudopotentials decay with in- 
creasing relative angular momentum. For spin and val- 
ley polarized electron systems this decay is determined 
by the ratios v/^V^^i"-* and V:^^"''' /vt\ A more stable 
FQHE is expected for Landau levels when the ratio of 
the pseudopotentials is large. 

In Table U the values of the ratios are shown for 
graphene and for non-relativistic systems for two low- 
est Landau levels with n = and 1 (only in these Lan- 
dau levels the FQHE can be observed). For the non- 
relativistic system the most stable FQHE is observed for 
the n = Landau level, which is supported by the data 
in Table U where ^/^^"■'/Vg''"' for the non-rclativistic sys- 
tem is the largest in the n — Landau level. A different 
situation occurs for the graphene system. Here for the 
n = Landau level, the wave functions are identical to 
the non-relativistic n = Landau level wave functions 
[Eqs. (HI)-®]. Therefore, the properties of the FQHE 
for the n — Landau levels of a non-relativistic sys- 
tem and graphene are the same if expressed in units of 
Efj. The wave functions of the n = 1 Landau level of 
graphene is the mixture of the n = and n ~ I non- 
relativistic wave functions, which results in an enhance- 
ment of the electron-electron interaction strength for the 
n = 1 graphene Landau level [13 • this Landau 
level, the ratio v}"'' /V^"^ has the largest value (Table HI, 
which suggests that the gaps of the FQHE states should 
have the largest value in graphene for the n — 1 Landau 
level. 

In theoretical studies the FQHE is often analyzed by 
numerical diagonalization of the Hamiltonian matrix for 
finite-size electron systems in either the planar (torus) or 
the spherical geometry [l9| . In the case of the spheri- 
cal geometry [3, [l^ [20] the magnetic field is introduced 
in terms of the integer number 25 of magnetic fluxes 
through the sphere in units of the flux quantum, where 
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= (graphene) 
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1.26 
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= 1 (graphene) 


1.68 


1.33 


n 


= 1 (non-relativistic) 


1.32 


1.36 



the radius of the sphere R is deflned as i? = \/SiQ . The 
number of available states in a sphere, which corresponds 
to the states of a single Landau level in planar geome- 
try, is (25* -I- 1). For a given number of electrons the 
parameter S determines the filling factor of the Landau 
level. For example, for the filling factor ;/ = 1/m it is 
S = (m/2)(7Vg — 1). In the case of the many-electron 
system the lowest eigenvalues of the interaction Hamil- 
tonian matrix determine the nature of the FQHE state 
and the corresponding neutral excitation gap [21|. The 
numerical results obtained for a finite size system show 
that the FQHE excitation gap i n gr aphene is the largest 
for the ri = 1 Landau level [13, l23 - l24| . As an example, 
for — 8 electrons the excitation gap is O.OSSe^^ for the 
n = Landau level and 0.094e^ for the n = 1 Landau 
level. This behavior is consistent with the properties of 
the pseudopotentials shown in Table H] 

In Fig. [T] the Landau levels in graphene correspond- 
ing to indices n — 1 and -1 and having the strongest 
electron-electron interactions, which results in a more 
stable FQHE, are shown by red lines. The electron- 
electron interaction in the Landau level with index n = 0, 
shown by a blue line, is identical to the interaction in the 
n = Landau level of the non-relativistic system. Exper- 
imental observation of FQHE in a suspended graphene 
[2B-I23 and robustness of the FQHE plateaus, which were 
observed even at a weak magnetic field ~ 2 Tesla, illus- 
trate the enhancement of the electron-electron interac- 
tions in graphene as compared to that for a conventional, 
non-relativistic semiconductor systems. 



III. INTERACTING FERMIONS IN BILAYER 
GRAPHENE 

Additional control of the interaction properties of 'rel- 
ativistic' (Dirac-like) particles in graphene is also possi- 
ble in a system of bilayer graphene [29l - l33j which con- 
sists of two coupled graphene layers. Bilayer graphene 
has been intensely investigated because of its intriguing 
properties. The effective low-energy Hamiltonian in this 
case is similar to the Dirac-like nature of that in mono- 
layer graphene, but with a quadratic (instead of linear) 
off-diagonal term [4, ^]. The low-energy dispersion is 
also quadratic. The massive Dirac fermions in bilayer 
graphene also posses the pseudospin degree of freedom 
and are chiral [3, [3^ . In a perpendicular magnetic field 
the Landau levels in a bilayer graphene follow the se- 
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FIG. 2: Schematic illustration of two different types of stack- 
ing of bilayer graphene, consisting of two coupled monolay- 
ers of graphene: (a) AA stacking; (b) Bernal stacking. Each 
graphene layer consists of two inequivalent sites A and B. The 
intra-layer and intra-layer hopping integrals are shown by 70 
and 7i , respectively. 



quence, £„ cx \/\n\ {\n\ — 1 ) fo r n > 1 with a doubly 
degenerate Eq = for n = pol . IsBj . An important char- 
acteristic of bilayer graphene is that it is a semiconductor 
with a tunable bandgap between the valence and conduc- 
tion bands [s^ . This property modifies the Landau level 
spectrum and influences the role of long-range Coulomb 
interactions fs^l • In a magnetic field, the electronic prop- 
erties of the graphene bilayer can be controlled by (i) the 
inter-layer bias voltage applied to two graphene monolay- 
ers, (ii) the intra-layer asymmetry potential due to the 
contact of one of the layers with a substrate, (iii) by ap- 
plying an in-plane magnetic field, and (iv) by introducing 
mechanical deformation of bilayer graphene [s^ . Below 
we consider only the effects of a bias potential and an 
in-plane magnetic field. 



IV. BIASED BILAYER GRAPHENE 

Bilayer graphene comprises two coupled graphene 
monolayers [30|. Depending on the orientation of the 
monolayers, there are two main stacking of a graphene 
bilayer: (i) the AA stacking and (ii) the Bernal (AB) 
stacking, which are shown schematically in Fig. [21 There 
is also the intermediate type of stacking of two mono- 
layers corresponding to the rotated bilayer graphene, in 
which monolayers are rotated relative to each other by 
an arbitrary angle [s^-H^ . These systems show rich low- 
energy physics due to the modulated nature of the inter- 
layer coupling. 

For the AA stacking [Fig. 2 (a)] in a perpendicular 
magnetic field, the interlayer coupling occurs between the 
Landau levels of the two layers with the same Landau 
level indices. This coupling changes the energies of the 
Landau levels of the monolayers, but does not affect the 
wave functions of the layers. Therefore, the pseudopoten- 
tials, which characterize the electron-electron interaction 
properties, are completely identical to the corresponding 
pseudoptentials of a monolayer graphene. In this case the 
FQHE in the Landau levels of bilayer graphene with AA 
stacking has the same properties, e.g., the same FQHE 



gaps, as in the corresponding graphene monolayers. 

In the case of Bernal (AB) stacking [Fig. 2 (b)], the 
interlayer coupling strongly modifies the properties of the 
Landau levels in the system. The Hamiltonian of the 
graphene bilayer with AB staking for valley ^ = ±1 can 
be written as [Sll 
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(11) 



where U is the inter-layer bias voltage and "fi « 
400 meV is the interlayer hopping integral. The 
corresponding wave function is described by a four- 
component spinor (V'Aj ' V'b^ i V'b^' ^Aij)"^ valley K 
and (V'b2''/'a2'V'Ai'V'b^)"^ for valley K' . Here the sub- 
indices A-^^, B^, and A2, B2 correspond to lower and upper 
layers respectively. The wave function corresponding to 
the Hamiltonian dTTI) has the form 



C2 



(12) 



where C^, Cj, C3, and C4 are constants. Therefore, the 
wave functions in bilayer graphene with Bernal stacking 
is a mixture of the non-relativistic Landau wave functions 
with indices n — 1, n, and n -\- 1. 

In the expression for the wave function (|12p of bilayer 
graphene, the Landau index n can take the following val- 
ues: n = —1,0,1,.... Here we assume that if the index 
of the Landau wave function is negative then the func- 
tion is identically zero, i.e., 4>-2m = ^^id = 0. 
In this case, for n — —1 the wave function is just 

^-''i!™ = (0,0,0, (/.„,,„), i.e., the coefficients C^, C^, C3 
are equal to zero. There is only one energy level corre- 
sponding to n = — 1. For n = 0, the wave function ([T2|) 
has a zero coefficient C'l, which results in three energy 
levels corresponding to n = 0. For other value of n, i.e., 
for n > 0, there are four eigenvalues of the Hamiltonian 
(fTTj) . corresponding to four Landau levels in a bilayer for 
a given valley ^ = ±1. The eigenvalue equation deter- 
mining these Landau levels, have the form 36] 

(£ + ^uf - 2nj \(e - ^uf - 2{n + 1)1 = 7^ [e^ - u^] , 

(13) 

where e is the energy of the Landau level in units of e ^ 
{eg = fivp/eg), u = U/{2eg), and 71 = 7i/es- ^ is 
convenient to introduce the following labeling scheme for 
the Landau levels determined by Eq. (|13p . The four Lan- 
dau levels correspond to two valence levels which have 
negative energies, and two conduction levels, which have 
positive energies. Then the four Landau levels of bilayer 
graphene for a given value of n and a given velley ^ can 
be labelled as n^j^\ where z = — 2,— 1,1, 2 is the label of 
the Landau level in the ascending order. Here negative 
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and positive values of i correspond to the valence and 
conduction levels, respectively. The Landau levels of dif- 
ferent valleys are related through the following equation 

^('^i^^) — Although for n = there are only 

three Landau levels and for n — —1 there is only one 
Landau level, it is convenient to include the n = — 1 Lan- 
dau level into the set of n = Landau levels and label 
them as of \ where i — —2, —1, 1, 3. 

At the zero bias voltage, the Landau levels become two- 
fold valley and two-fold spin degenerate and are given by 
the expression 



±\ 2n- 



1 



l + -^±-^(2 + 7?)2 + 8n7f. (14) 



Since the FQHE is expected only in the Landau levels 
with low values of the index, n, we consider below the 
sets of Landau levels of bilayer graphene with ?i = and 
n = 1 only. The wave functions of these Landau levels 
are mixtures of the conventional, non-relativistic Landau 
functions with indices 0, 1, and 2. 

Once the wave functions (fT^ of the bilayer Landau lev- 
els are obtained, the form factor in the pseudopotentials 
(III) can be obtained from 

F„(g) = |CipL„_i(gV2)+(|C2p + |C3p)L„(gV2) 
+ \C^\^L^^^[qV2). (15) 

With the known form factors, the pseudopotentials, 
which determine the interaction strength and the FQHE 
in a given Landau level, can be calculated. 

There are two special Landau levels of bilayer 
graphene. For n ~~\ there are two solutions (one for 
the valley K and one for K') of Eq. ([T3| with energies 
£ = — ^w. The corresponding wave function 



(bi) 



(bi) 









(16) 



is determined only by the n = non-relativistic Landau 
level wave function. Therefore the FQHE and the inter- 
action properties of these Landau levels are exactly the 
same as those for the 0-th conventional (non-relativistic) 
Landau level. 

For n = and for small values of u there is another 
solution of Eq. (|13p with almost zero energy, e « 0. The 
wave function of this Landau level has the form 



m ' ,m 



( r \ 









(17) 



For a small magnetic field, Ji, the wave function 

becomes (V'l 0, 0, 0)-^ and the Landau level becomes 
identical to the n = 1 non-relativistic Landau level. In a 
large magnetic field £^ 7^ , the Landau level wave func- 
tion becomes (0, 0, ipg O)"'" and the bilayer Landau level 
has the same properties as for the n = non-relativistic 
Landau level. The corresponding form factor of the Lan- 
dau level (IT7)) is given by 



^0 , (9) 



7? 



Li(gV2) 



2f2 



24 



io(9V2). 



\ 7l01,m / 



(18) 

With increasing magnetic field, i.e., with increasing e^, 
the bilayer Landau level Q_i becomes identical to (i) the 
n — \ non-relativistic Landau level with the form factor 
of Li{q^ /2) for small B, <C 7i; (ii) the n = 1 Landau 
level of the monolayer graphene with the form factor of 
i[Lo((zV2) + Li((7V2)] for = 7i/V2; and, (iii) the 
n — Q non-relativistic Landau level with the form factor 
of L^lq-^ /2) for large B, Sg 7^. For typical values 
of the interlayer coupling, — 400 meV, the condition 
= 72 7-^/2 is achieved for a large magnetic field B = 
120 Tesla. Under this condition only the first regime can 
be experimentally realized. Below we show that an in- 
plane magnetic field can suppress the interlayer coupling, 
which opens the possibility of experimental observation 
of transitions between the different regimes (i)-(iii). 

For all the Landau levels [except the levels described by 
Eq. (jl6l) ] in bilayer graphene the electron-electron inter- 
action strength and the stability, i.e., the excitation gaps 
of the FQHE states depend on the magnetic field B, the 
bias voltage U, and the Landau level index. Therefore 
the interaction properties of a bilayer graphene can be 
controlled by the external parameters (43| , which is a dif- 
ferent situation than in a monolayer graphene, where the 
interaction properties depend only on the Landau level 
index. 

The stable FQHE states in a bilayer graphene are ex- 
pected for the n = and n — 1 Landau level sets, which 
are the mixtures of the n = 0, n = 1, and n = 2 non- 
relativistic Landau level wave functions. This mixture 
depends on the values of the parameters of the system. 
To characterize the stability of the FQHE we evaluate 
numerically the FQHE excitation gaps for a finite size 
system in a spherical geometry. We present below the 
results for the v = ^ FQHE state. A similar behavior 
is expected for other main fractions of the FQHE, i.e., 

~ 5' 5 ' 3 

In Fig. [3] we show the dependence of the Landau levels 
on the bias voltage U for a fixed magnetic field and for 
different valleys. The results are presented only for the 
Landau levels with indices n — Q and 1, i.e. only for the 
Landau level where the FQHE can be observed. The cor- 
responding — \ FQHE gaps are shown in Fig. [3] (b,d). 
Both for the K and K' valleys there is a special Landau 
level, 0^+^ (for the K valley) and O^Zi (for the K' val- 
ley) , that is described by the wave function of the type as 
in Eq. In these Landau levels the FQHE gap does 
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FIG. 3: The Landau levels of the bilayer graphene [panels (a) 
and (c)] shown as a function of the bias voltage, U. Panels (b) 
and (b): the Coulomb gaps of the FQHE in corresponding 
Landau levels. The results are obtained in spherical geometry 
for a finite-size system with eight electrons and 25* — 21 flux 
quanta. The numbers next to the lines are the labels of the 
Landau levels. The same type of lines [in panels (a) and 
(b) and panels (c) and (d)] correspond to the same Landau 
levels. Panels (a) and (b) correspond to the valley K' , while 
panels (c) and (d) correspond to the valley K. The system is 
characterized by 7^ — 400 meV and a magnetic field _B = 15 
Tesla. The arrows in panels (b) and (d) indicate the gap of the 
i— FQHE in the n = 1 Landau level of a monolayer graphene. 



not depend on the bias voltage and is exactly equal to the 
FQHE gap of a conventional (non-relativistic) semicon- 
ductor systems for the n = Landau level. In all other 
levels the FQHE gap depends on the bias voltage, which 
clearly illustrates the sensitivity of the interaction prop- 
erties on the external parameters, i.e., the bias voltage. 
Although the interaction strength within a single Landau 
level can be controlled by the bias voltage, the results il- 
lustrated in Fig. [3] show that the FQHE gaps in bilayer 
graphene are usually less than the largest FQHE gap in 
a monolayer graphene. This FQHE gap in a monolayer 
graphene is realized in the n — 1 Landau level and is 
shown by red arrows in Fig. [3l 

For a smaller interlayer hopping integral, the Landau 
levels in bilayer graphene show anticrossings as a func- 
tion of the bias voltage [isj. These anticrossings result 
in a strong mixture of different Landau levels, which can 
greatly modify the properties of the Landau level wave 
functions and change the interaction strength within a 
single Landau level. This behavior is illustrated in Fig. 
m where the dependence of the Landau levels on the 
bias voltage is shown for = 30 meV. The anticross- 
ings of the Landau levels result in transitions from an 
incompressible state (FQHE) to a compressible state (no 
FQHE) within 

a single Landau level (see the Landau level in Fig. 
13] (a)). There is also a double transition, marked by the 
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FIG. 4: A few lowest Landau levels of the conduction band 
as a function of the bias potential, U, for inter-layer cou- 
pling of 7i — 30 meV and a magnetic field of 15 Tesla. The 
numbers next to the curves denote the corresponding Landau 
levels. Panels (a) and (b) correspond to the K and K' val- 
leys, respectively. The Landau levels where the FQHE can 
be observed are drawn as blue and green filled dots. The 
green dots correspond to the Landau levels where the FQHE 
states are identical to that of a monolayer of graphene or a 
non-relativistic conventional system. The red dots represent 
Landau levels with a weak FQHE. The open dots correspond 
to Landau levels where the FQHE is absent. In (a), the dashed 
lines labeled by (i) illustrates the transition between FQHE 
(symbol 'F') and no FQHE (symbol 'NF') states under a con- 
stant gate voltage and variable bias potential [i^ . 



dashed line (i), at the Landau level l^'''''. At this Landau 
level, the electron system with fractional filling shows 
transitions from an incompressible state (FQHE) at small 
bias voltage J7 to a compressible state (no FQHE) at 
intermediate values of U and then to an incompressible 
state (FQHE) at large U. No such transition has ever 
occured in conventional semiconductor systems. 

Although for experimentally realized bilayer systems 
the interlayer hopping integral is relatively large, ~ 
400 meV, the interlayer coupling can be controlled and 
suppressed by an applied in-plane magnetic field. This 
situation is discussed in the next section. 



BILAYER GRAPHENE IN A TILTED 
MAGNETIC FIELD 



A tilted magnetic field, applied to a quasi-two- 
dimensional electron system, can modify the electron dy- 
namics and correspondingly the electron-electron interac- 
tion strength. In a graphene monolayer, due to its purely 
2D nature, the component of the magnetic field parallel 
to the monolayer does not influence the electron's spatial 
dynamics, although it can alter the electron spin dynam- 
ics, which is sensitive to the total magnetic field p4. |45|. 

Bilayer graphene is a quasi-two-dimensional system. 
The electron dynamics in such a system is sensitive to 
both perpendicular and in-plane components of the mag- 
netic field [46]. To introduce a tilted magnetic field 
into the Hamiltonian of bilayer graphene, we introduce 
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the vector potential A = {0, Bj_x, B^^y), where Bj^ {z- 
component) and B^^ (x-component) are perpendicular 
and in-plane components of the tilted magnetic field. 
Here the z axis is perpendicular to the graphene monolay- 
ers. The perpendicular component of the magnetic field, 
which alters the electron dynamics in the (x, y) plane, 
is introduced in the bilayer Hamiltonian by replacing 
the 2D momentum by the generalized momentum [simi- 
lar to the Hamitonian ([TT])]. The parallel component of 
the magnetic field is introduced through the Peierls sub- 
stitution by multiplying the interlayer hopping integral 
7^ with the phase factor exp{—ie/hA^d) = exp(— i/3y), 
where /3 = eB^^d/h, and d is the interlayer separation. 
Then the Hamiltonian of the bilayer graphene with AB 
stacking and at zero bias voltage becomes 



n 



{AB) 



1 





Wp7r_ 





' \ 


































Wp7r+ 


1 



(19) 



The wave functions corresponding to the Hamiltonian 
([T9| can be expressed in terms of the non-relativistic 
2D Landau wave functions. For the vector potential 
Aj^ = (0, i?j^x,0) corresponding to the perpendicular 
component of the magnetic field, the 2D Landau wave 
functions are parametrized by the y component of the 
wave vector and the Landau index n, and are described 
as 



X exp 



211 



Here the 



where V/c(a;) = C'o^ini (t^") "xp 
are the Hermite polynomials, and x^. 
magnetic length f q is defined by the perpendicular com- 
ponent of the magnetic field, ~ ^hjeB^. Then the 
wave functions of the Hamiltonian (jl9p are parametrized 
by the Landau level index n and the wave vector fc and 
have the form 



/ ^C'i'^|„|-i,fc \ 



(bi: 



^|n|,fc-/3 
3<^|n|,fe-|-/3 
V CC4</'|„| + i,fe } 



(21) 



where Cj^, C^-, C'g, and C4 are constants. The in-plane 
component of the magnetic field results in a coupling of 
the Landau wave functions with the wave vectors fc, fc— /3, 
and fe -I- /3. For the wave functions of the form (l?T|) the 
Hamiltonian of bilayer graphene in a Landau level with 
index n takes the form 



(AB) 








Up7r_ 
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(22) 
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FIG. 5: The Landau levels of bilayer graphene in a tilted 
magnetic field and zero bias voltage shown as a function of the 
parallel component of the magnetic field. The perpendicular 
component of the tilted magnetic field is (a) 1 Tesla and (b) 2 
Tesla. The labels next to the lines denote the corresponding 
Landau levels, where only the Landau levels in which the 
FQHE can be observed, are labelled. Only the Landau levels 
with positive energies are shown. 



where 



'«„(m) = / dxtPf^{x)ij^^_i^{x) 



(23) 



depends on the dimensionless parameter /i = I3£q ~ 
eB^^diQ/h. Therefore the effect of the in-plane component 
of the magnetic field on the electron dynamics in a bi- 
layer graphene is the reduction of the interlayer coupling, 
i.e., the interlayer coupling is 7^^ „ = 7i'*n(M) < 7i- This 
reduction depends on the Landau level index and on the 
dimensionless parameter fi. Due to the small interlayer 
distance d — 3.3 A, the parameter ^ is relatively small. 
To increase the value of this parameter, the perpendic- 
ular component of the magnetic field needs to be small, 
i.e., the magnetic length should be large, and the paral- 
lel component of the magnetic field should also be large. 
For example, for Bj^^l Tesla, n = 0.014Bj| [Tesla]. For 

the first lowest Landau level indices the function K^^^fJ-) 
is Ko(m) = e-^'/4 and n^ifi) ^ e-^'/* (l - 

The Landau levels of a bilayer graphene in a tilted 
magnetic field are given by Eq. (jl4[) . in which the inter- 
layer coupling 7]^ should be replaced by 7]^ „ — n/^B- 
Here is calculated in terms of the perpendicular com- 
ponent of the tilted magnetic field. In Fig. [S]the depen- 
dence of the Landau levels on the parallel component of 
the magnetic field, , is shown for a few lowest Lan- 
dau levels of a bilayer graphene. Increasing the paral- 
lel component of the magnetic field, the energies of the 
Landau levels are reduced, which is consistent with the 
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VI. THE PFAFFIANS IN GRAPHENE 



to 0.084 
O 
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FIG. 6: The FQHE gaps at different Landau levels of bilayer 
graphene shown as a function of the parallel component of the 
tilted magnetic field. The perpendicular component of the 
magnetic field is 1 Tesla. The labels next to the lines denote 
the corresponding Landau levels [Fig. [5]. The bias voltage 
is zero in this case. The results are obtained in spherical 
geometry for a finite-size system with eight electrons and 2S = 
21 flux quanta. 



For conventional 2D semiconductor systems, in ad- 
dition to the usual incompressible FQHE states that 
are realized for the odd-denominator filling factors, the 
electron-electron interaction is also responsible for the 
formation of a special type of incompressible state cor- 
responding to the fractional filling factor = §■ This 
filling factor corresponds to a completely occupied n = 
Landau level with two spin components and a half-filled 
n = 1 Landau level. Since the completely occupied Lan- 
dau levels do not contribute to the electron dynamics, 
the ground state of the ^ = \ system is determined 
by the electron-electron interaction alone in the n — 1 
Landau level. The incompressible state with a large ex- 
citation gap is formed in this half-filled Landau level. 
One unusual property of this state is that the elementary 
charged excitations have a charge e* = e/4. They obey 
the "non-abelian" statistic s |4 8L l49l | and carry the signa- 
ture of Major ana fermions [50[. It was proposed that the 
ground state of the half-filled n = 1 Landau level is de- 
scribed by a Pfaffian j^^j HH or the anti-Pfaffian function 
[52,^1. The Pfaffian state is written 



reduction of the interlayer coupling, ^-^ with increas- 
ing _B||. The dependence of the Landau levels on By 
becomes weaker with increasing perpendicular magnetic 
field [Fig. [5] (a,b)]. Therefore, the effect of an in-plane 
magnetic field on the Landau levels can be observed only 
for a small perpendicular magnetic field, Bj^ « I, and a 
large parallel magnetic field, > 50 Tcsla. It should be 
pointed out that although the perpendicular component 
of the field is rather small, in a conventional semiconduc- 
tor system the FQHE has been reported in a magnetic 
field of S < 3 Tesla 113]. 

The interaction properties of electrons in the Landau 
levels of a bilayer graphene also depend on the in-plane 
component of the magnetic field. This dependence is 
visible only for small perpendicular components of the 
magnetic field, i.e., Bj^ « I Tesla. In Fig. [5] we show the 
-|-FQHE gap as a function of the in-plane component of 
the magnetic field for different Landau levels. For Bj^ « 1 
Tesla, only three Landau levels (with positive energies) 

support the FQHE states. One Landau level O^^', the 
wave function of which has the form of (|16p and depends 
only on the perpendicular component of the magnetic 
field, does not show any dependence on the in-plane com- 
ponent of the magnetic field. The interaction strength in 
the Landau levels and depends weakly on B^^ 
[Fig. [S]. The interaction strength increases with B^^ for 

the Landau level 02^'' and decreases for the Landau level 

l^^''. Therefore the parallel component of the magnetic 
field can in fact, enhance the electron-electron interac- 
tion strength for some Landau levels (Fig. [5]) in a bilayer 
graphene. 



^-pf - Pf I 



V i<j 



(24) 

where the positions of the electrons are described in terms 
of the complex variable z — x ~ iy and the Pfaffian is 
defined as [2a, LI] 



Pf Af, = 



1 



N/2 



y 2N/2 (7V/2)! 



J2 sgnaHK 



(2/-1)ct(2;) I 



•res. 



1=1 



(25) 

for an N X N antisymmetric matrix whose elements are 
Mjj. Here Sj^ is the group of permutations of N objects. 

The Pfaffian state is the exact ground state with zero 
energy for the electron system at half filling with a spe- 
cial three-particle interaction which is non-zero only if all 
three particles are in close proximity to each other [s^ j. 
For realistic two-particle interactions the Pfaffian state 
is not an exact eigenstate of the half-filled system. In 
the case of the Coulomb interaction, the overlap of the 
ground state of the v — ^ system in the n = 1 Landau 
level with the Pfaffian function is around 80%. By vary- 
ing the two-particle interaction potential, i.e., the pseu- 
dopotentials, a stronger overlap (~ 99%) of the ground 
state of the v — \ system in the n — 1 Landau level with 
the Pfaffian state is possible. The proximity of the actual 
V — \ ground state to the Pfaffian state is most sensi- 
tive to the lowest pseudopotentials, V^, V^, and Vr^. In 
graphene, there are two lowest Landau levels with indices 
n — Q and n = 1 with strong electron-electron interac- 
tions. Although conventional FQHE states with odd- 
denominator filling factors can be observed at these Lan- 
dau levels, the Pfaffian state with half-filling of the cor- 
responding Landau level cannot be realized (ssj . In the 
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n = Landau level in graphene, the interaction potential 
is identical to the one in the n — non-relativistic Lan- 
dau level, and similar to the case of the non-relativistic 
system, the PfafBan state is not the ground state of the 
= i system in the n = Landau level in graphene. 
In the n = 1 Landau level in graphene, although the 
electron-electron interaction results in the stable odd de- 
nominator FQHE states, the ground state of the half- 
filled Landau level is compressible and is not described 
by the Pfaffian function. The overlap of the ground state 
of the 1^ — system with the Pfaffian function is less than 
0.5 in all Landau levels of the monolayer graphene [s^ . In 
bilayer graphene the interaction strength and the corre- 
sponding Haldane pseudopotentials can be controlled by 
the external parameters, such as the bias voltage and the 
direction of the magnetic field. In this case the stability 
(i.e., the magnitude of the excitation gap) of the Pfaf- 
fian state can be strongly enhanced. In bilayer graphene 
there are two 'special' Landau levels (for valley K) 

and oj"'' (for valley K'), that are described by Eq. (IT71) . 
The numerical calculations in a spherical geometry show 
that only at these special Landau levels the overlap of 
the ground state with the Pfaffian state is large [SS?]. In 
all other bilayer Landau levels the overlap of the = ^ 
ground state with the Pfaffian state is found to be small 
(< 0.6) and these states cannot be described by the Pfaf- 
fian function. In spherical geometry, the Pfaffian state 
in a system of electrons is realized for the parameter 
25' = 2N^ — 3, which corresponds to filling factor = \ 
in the thermodynamic limit. 

For the zero bias voltage, the energies of the Landau 
levels 0L"^i^ and 0^"' [Eq. ^] are zero and the levels 
are degenerate with zero-energy Landau level given by 
Eq. (|T6)) . For a finite bias voltage, the degeneracy of the 
Landau levels is lifted. The form factor in the Lan- 
dau level 0^^ is calculated from Eq. and determines 
the interaction properties of the electron system in the 
level 0^''. For a small magnetic field, ^ e^, the form 
factor is identical to the form factor of the non-relativistic 
n — 1 Landau level. Therefore, in this limit the ground 
state of the v = half-filled system in the O^"* Landau 
level is incompressible and is determined by the Pfaffian 
state. In a large magnetic field, 7^ eg, the form factor 
becomes identical to that of the n = non-relativistic 
system, for which the = 5 state is compressible. For in- 
termediate values of the magnetic field, the — system 

in the 0^' Landau level shows an unique behavior as a 
function of the magnetic field: with increasing magnetic 
field the overlap of the ground state of the system with 
the Pfaffian state shows a maximum for a finite value of 
the magnetic field [55j. Therefore, the stability of the 
Pfaffian = ^ state in bilayer graphene can be increased 
when compared to that in non-relativistic systems. 

Our results shown in Fig. [7] illustrate the non- 
monotonic dependence of the interaction properties of the 
V — \ system in the ol!*^-* Landau level. Here the over- 



400 meV 




10 15 20 25 30 
Magnetic Field (T) 



FIG. 7: (a) Overlap of the exact many-particle ground state 
with the Pfaffian function, (b) Collective excitation gap of 
the I' ~ \ state. The results are for A'^^ = 14, 2S — 25, and 
[7 = 5 meV. The black and red lines correspond to 71 = 400 
meV and 300 meV, respectively. The results are shown for 



the = I system in the Landau level 0^'. 



lap of the ground state with the Pfaffian state and the 
corresponding excitation gap are shown. With increasing 
magnetic field the properties of the system change non- 
monotonically and for 7^^ — 400 meV the overlap with the 
Pfaffian state has a maximum in a magnetic field of ~ 10 
Tesla. The corresponding excitation gap also has a max- 
imum at this point. In dimensionless units the maximum 
appears when ^i/e^ « 4.9. Therefore, for smaller values 
of "fi , the maximum of the overlap is realized for a smaller 
value of magnetic fields (see the results for -fi — 300 meV 
in Fig.IZl). 

It is possible to suppress the interlayer hopping integral 
by applying a tilted magnetic field where the in-plane 
component of the magnetic field determines the suppres- 
sion of 7]^ by a factor of k„ . To identify the effect of the 
in-plane magnetic field on the stability of the Pfaffian 
state, we characterize the interaction properties of the 
half-filled Landau level by the ratios of the pseudopoten- 
tials corresponding to the lowest relative angular momen- 
tum. The V = ^ Pfaffian state is most sensitive to two 
parameters of the pseudopotentials V^/Vg and V^a/Vg [56j . 
In Ref. [ssj, in the thermodynamics limit, the region of 
the pseudopotenial parameters, for which the most sta- 
ble Pfaffian state can be realized, was obtained in the 
plane (V^/Vg) — iV^/Vr^). We apply that approach on 
our bilayer system and evaluate the pseudopotential pa- 
rameters Vi/V^ and V^/V^ in the Landau level 0^^ as a 
function of the magnetic field. In this way, we can iden- 
tify the regions of the magnetic field with the most stable 
Pfaffian state. In Fig. |S] the parameters Vi/Vc^ and 
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FIG. 8: Ratios of pseudopotentials for two values of the an- 
gular momentum V1/V5 [panel (a)] and 14/14 [panel (b)] as 
a function of the perpendicular component of the magnetic 
field and for different parallel components of the magnetic 
field, S|| = 5, 50, and 100 Tesla. The data are shown for the 

Landau level 0^' . The hatched regions correspond to the val- 
ues of the pseudopotentials where one observs a large overlap 
of the ground state with the Pfaffian state and also a large 
excitation gap of the incompressible ground state. 

are shown as a function of the perpendicular component 
of the magnetic field and for different in-plane compo- 
nents of the magnetic field. These results demonstrate 
that with increasing parallel component of the magnetic 
field the values of the pseudopotentials, which correspond 
to the most stable Pfaffian state and which are illustrated 
by the hatched region in Fig. [H are obtained for smaller 
values of the parallel magnetic field. Therefore, for a 
given value of the perpendicular component of the mag- 
netic field the parallel component of the magnetic field 
increases the stability of the Pfaffian state. 

Another interesting effect, introduced by the parallel 
component of the magnetic field, is the strong modifica- 
tion of the interaction properties of the electron system 
at small values of the perpendicular magnetic field. Such 
changes in the interaction potential result in an enhance- 
ment of the stability of the Pfaffian states for small values 
of As an example, for By — 100 Tesla the stability 
of the Pfaffian state is strongly increased for a weak per- 
pendicular magnetic field, < 2 Tesla (Fig. [5]). 



VII. CONCLUDING REMARKS 

In a magnetic field the strength of the electron-electron 
interaction, which is characterized by the value of the 
FQHE gap in a given Landau level, depends on the 
Landau level index and the external parameters of the 
graphene system. In the case of monolayer graphene, 
there are two types of Landau levels with indices n = 
and n = 1, which have strong electron-electron interac- 
tions, i.e., the FQHE can be observed only in these Lan- 
dau levels. Among these Landau levels, the strongest 
electron-electron interactions are realized in the n = 1 
graphene Landau level. In the Landau level with index 
n = 0, the interaction strength is exactly the same as 
that in the n = Landau level of the conventional (non- 
relativistic) system, which results in exactly the same 
FQHE gaps. 

The strength of the electron-electron interaction can 
be further controlled in a bilayer graphene, where the ad- 
ditional parameters that govern the interaction strength 
are the interlayer coupling, the bias voltage, and the ori- 
entation of the magnetic field. The bias voltage between 
the graphene monolayers changes the structure of the 
wave functions of the Landau levels and can strongly 
modify the interaction strength. In some Landau lev- 
els the electron-electron interaction can be stronger than 
that in a monolayer graphene, resulting in a more sta- 
ble FQHE. In a given Landau level and as a function of 
the bias voltage the bilayer graphene system can show 
transitions from a state with weak electron-electron in- 
teraction (FQHE being absent) to a state with strong 
electron-electron interactions (presence of FQHE). In a 
bilayer graphene, the electron-electron interactions can 
be additionally controlled by the direction of the mag- 
netic field, i.e., in a tilted magnetic field. The sensitivity 
of the interaction strength to the parallel component of 
the magnetic field is visible only for a weak perpendicular 
component of the magnetic field, sa 1 Tesla and for a 
strong parallel component of the magnetic field, > 50 
Tesla. Finally, we describe the stability of the Pfaffian 
state and the excitation gap in a half-filled n—\ Landau 
level in bilayer graphene. We also discuss the possibility 
of making the Pfaffian state more stable by applying a 
tilted magnetic field. 
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